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New data are reported from the operation of a 4.0 kg CF3I bubble chamber in the 6800-foot- 
deep SNOLAB underground laboratory. The effectiveness of uhrasound analysis in discriminating 
alpha-decay background events from single nuclear recoils has been confirmed, with a lower bound 
of >99.3% rejection of alpha-decay events. Twenty single nuclear recoil event candidates and three 
multiple bubble events were observed during a total exposure of 553 kg-days distributed over three 
different bubble nucleation thresholds. The effective exposure for single bubble recoil-like events was 
437.4 kg-days. A neutron background internal to the apparatus, of known origin, is estimated to 
account for five single nuclear recoil events and is consistent with the observed rate of multiple bubble 
events. The remaining excess of single bubble events exhibits characteristics indicating the presence 
of an additional background. These data provide new direct detection constraints on WIMP-proton 
spin-dependent scattering for WIMP masses >20 GeV/c^ and demonstrate significant sensitivity 
for spin-independent interactions. 

PACS numbers: 29.40.-n, 95.35. -|-d, 95.30.Cq 



I. INTRODUCTION 

There is abundant evidence that '^85% of the matter 
in the Universe is cold, dark, and nonbaryonic The 
leading candidate for the dark matter is a relic density, 
left over from the big bang, of an as yet undiscovered 
weakly interacting massive particle 0- If weakly inter- 
acting massive particles (WIMPs) are the dark matter, 
then they may scatter off nuclei with enough energy and 
at a high enough rate to be detectable in the laboratory 
through the observation of single recoiling nuclei [31] . 

The Chicagoland Observatory for Underground Par- 
ticle Physics (COUPP) employs a novel bubble cham- 
ber technique to search for the single nuclear recoils 
that would arise from WIMP-nucleon elastic scattering 
The physics of bubble nucleation provides a pow- 
erful natural discrimination between nuclear recoils and 
the electron recoils from the abundant gamma-ray and 
beta-decay backgrounds. If the chamber pressure and 
temperature are chosen appropriately, electron recoils do 
not nucleate bubbles Nuclear recoil backgrounds in 
COUPP can still arise from neutron interactions or from 
the alpha decay of contaminants in the bubble chamber 
fluid. The chamber is surrounded by a low-Z water and 
polyethylene shield which moderates neutrons from spon- 
taneous fission and (alpha, n) in materials at the experi- 
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mental site to a negligible level. The 6800-foot (6010 m 
water equivalent) overburden of the SNOLAB site elim- 
inates neutrons of cosmogenic origin. Neutrons arising 
from detector materials interior to the shielding [6] can 
provide a limiting background, as discussed below. 

Because the bubble chamber is a threshold device with 
no event-by-event energy measurement, nuclear recoil 
events initiated by alpha decays provide a serious back- 
ground for a dark matter search. The use of acoustic dis- 
crimination has proven effective in mitigating the alpha- 
decay background 7, 8]. 

We report results from a 4.0 kg CF3I bubble chamber 
operated from September 2010 to August 2011 in the J- 
Drift of the SNOLAB deep underground laboratory. 
Results from the same bubble chamber, operated with a 
3.5-kg CF3I target in the MINOS underground area at 
Fermilab [Ifll were previously reported . 



II. EXPERIMENTAL METHOD 

The bubble chamber consisted of a 150-mm-diameter 
3-1 synthetic fused silica [ll| bell jar sealed to a flexible 
stainless steel bellows and immersed in propylene glycol 
within a stainless steel pressure vessel. The propylene 
glycol, which served as the hydraulic fluid to manage 
the inner pressure of the bubble chamber, was driven by 
an external pressure control unit. The flexible bellows 
served to ensure that the contents of the bell jar were 
at the same pressure as the hydraulic fluid, reducing the 
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stress in the silica vessel. The bell jar contained 4.0 kg 
of CF3I topped with water which isolated the CF3I from 
contact with any stainless steel surfaces or seals. The 
superheated CF3I was in contact only with the smooth 
synthetic silica surfaces or with the water interface above. 

The thermodynamic conditions of the chamber were 
monitored with two temperature sensors mounted on the 
bellows flanges and by pressure transducers which sepa- 
rately monitored the pressure of the hydraulic fluid and 
the inner vessel fluid. An additional fast ac-coupled pres- 
sure transducer monitored the pressure rise in the cham- 
ber to track bubble growth. Four lead zirconate (PZT) 
piezoelectric acoustic transducers epoxied to the exterior 
of the bell jar recorded the acoustic emissions from bub- 
ble nucleations, the audible "plink" used to trigger the 
flash lamps in early bubble chambers [l^l- Two video 
graphics array (VGA) resolution CCD cameras were used 
to photograph the chamber with a 20° stereo angle at a 
rate of 100 frames per sec. Stereo image data from the 
cameras were used to reconstruct the spatial coordinates 
of each bubble within the chamber. 

Each operating cycle of the bubble chamber began with 
the CF3I in its normal state, compressed to 215 psia. 
An expansion to the superheated state was accomplished 
by reducing the pressure from 215 psia to the operating 
pressure of 30.1 psia over a period of five sec. Following 
expansion and a 30-sec period for pressure stabilization, 
the chamber was live for the accumulation of dark matter 
data. In the expanded state, frame-to-frame differences 
in the image data provided the primary trigger for the ex- 
periment, typically initiating compression and capture of 
event data within 20 msec of a bubble nucleation. Com- 
pression and data capture were also initiated if consecu- 
tive pressure measurements indicated a possible bubble 
nucleation, if the operating pressure drifted out of the al- 
lowed range, if an error condition was detected, or if the 
chamber remained expanded beyond the 500-sec expan- 
sion timeout without a bubble nucleation. Return of the 
CF3I to its normal state under 215 psia compression was 
accomplished in 80 msec. The compression duration was 
30 sec, with a longer compression of 300 sec after every 
10th event to ensure that all CF3I gas produced by the 
bubble was condensed and returned to the liquid volume. 
During the compression period, the event data from the 
cycle were logged and the chamber was prepared for the 
next expansion. The expansion/compression cycle of the 
bubble chamber is illustrated in Fig. [T] Including the 57- 
sec average compression time, the 30 sec settling time, 
and the 500 sec maximum expansion time, the live-time 
fraction for the experiment could not exceed 84.4%. In 
practice the average live-time fraction ranged from 78.8% 
to 82.2% depending on the operating temperature of the 
chamber. 

The chamber was operated at a pressure of 30.1 psia to 
ensure good performance of the acoustic measurements. 
The bubble nucleation threshold was determined by the 
operating temperature of the CF3I. Dark matter search 
data were accumulated in three contiguous data sets at 
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FIG. 1: (Color Online) Pressure history from a sample event 
at 33.5° C. Time scale is linear within each region. The event is 
divided into four regions: (a) Chamber expands to the super- 
heated state, (b) Pressure regulation turns on at elapsed time 
of 5 sec and the chamber stabilizes by elapsed time 30 sec, 

(c) Chamber is live (accumulating dark matter data) from 
30 sec until a trigger or timeout at elapsed time of 500 sec, 

(d) Chamber compresses and sits compressed for 30 sec be- 
tween events, or 300 sec every tenth event. The mean ex- 
pansion times at 39.0°C, 36.2°C, and 33.5°C are 326, 396, 
and 417 sec, respectively. The shorter mean times at higher 
temperatures are due to an increased trigger rate during the 
expansion and stabilization periods. The majority of events 
at all temperatures end with a timeout. 

temperatures of 39.0°C, 36.2°C, and 33.5°C, correspond- 
ing to nominal Seitz model bubble nucleation thresholds 
of 7.8, 11.0, and 15.5 keV nuclear recoil energy [1^, re- 
spectively. To monitor the stability of the chamber, 12 
calibration runs with neutron sources were performed at 
scheduled intervals. Throughout the data taking, the 
performance of the chamber was stable and consistent 
with previous experience except for a higher rate of radon 
ingress into the active volume of the experiment. The 
rate of radon entering the active volume was '^8 atoms 
per day, resulting in 22 alpha-decay events per day, con- 
stant over the duration of the experiment. 

III. DATA ANALYSIS 

The reduction of the data consisted of examination of 
the photographic images to determine the number and 
spatial coordinates of bubbles, inspection of the pressure 
rise to confirm the bubble count and identify events oc- 
curring near the vessel walls, and analysis of the acous- 
tic traces to characterize the event types. Bubbles in 
the photographic images appear in sharp contrast to a 
retroreflective background and are identified in the image 
analysis algorithm as clusters of pixels that have changed 
significantly between consecutive frames. Reconstruction 
of the data from two stereo views provided the spatial co- 
ordinates of the bubble to a typical accuracy of a few mm, 
depending on the proximity of the bubble to the cameras. 
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The pressure rise analysis was based on data from an 
ac-coupled fast pressure transducer pl| which was sam- 
pled at 10 kHz for 160 msec around the onset of a nucle- 
ation. Empirically, the rate of pressure rise was well fit by 
a simple quadratic time dependence for bubbles formed 
in the bulk of the target fluid. The quadratic coefficient 
of the fit was found to be proportional to the number 
of bubbles in the event, and the quality of the fit was 
uniform over the volume of the experiment except near 
the boundaries. Because bubble growth is affected by 
the proximity of the bubble to the quartz vessel walls or 
the CF3I water interface, the quality of the quadratic fit 
deteriorated rapidly for bubbles near a boundary. The 
sensitivity of the bubble growth to the proximity of a 
boundary was studied using calibration neutron events 
where it was found that a modest cut on the chi-square 
reliably identified events that were near the vessel walls or 
the CF3I water interface. The pressure growth chi-square 
cut effectively provided a fiducial volume definition that 
was uniform around the perimeter of the chamber and 
performed somewhat better for this purpose than the 
stereo reconstruction of the camera images. The pres- 
sure growth fit was therefore used to provide the formal 
fiducial volume cut for the experiment. 

The third and final element of event reconstruction was 
the evaluation of the acoustic signals and classification 
of event types. The acoustic transducer signals were dig- 
itized with a 2.5-MHz sampling rate and recorded for 
40 msec for each event. The signals were filtered using a 
single-pole high-pass filter with a cutoff at 500 Hz, and a 
low-pass anti-aliasing filter cutting off at 600 kHz. The 
preevent baseline for each of the acoustic signals was ex- 
amined to determine the time of bubble formation, io- 
A fast Fourier transform was constructed for the times 
to — 1 msec < t < to + 9 msec. The sound of bubble nu- 
cleation showed a broad emission distinctly above back- 
ground noise up to a frequency of 250 kHz. The acoustic 
signature for a single recoiling nucleus was calibrated by 
studying events initiated by neutron sources. The acous- 
tic power was observed to vary slightly with the position 
of the bubble within the chamber, and the position de- 
pendence was found to vary with frequency. To account 
for the position and frequency dependence, the acoustic 
signal was analyzed separately in four frequency bands 
(1.5-12, 12-35, 35-150, or 150-250 kHz) which were sep- 
arately corrected for spatial dependence and normalized. 
The acoustic event discrimination was based on a single 
acoustic parameter AP which is a frequency weighted 
acoustic power density integral, corrected for sensor gain 
and bubble position. 



AP = AiT) J2 Gj ^ C„(x) £ / X psd'f. 



(1) 



where A{T) is an overall temperature dependent scale 
factor, Gj is the gain of acoustic transducer j, C„(a;) is 
the correction factor for the bubble position dependence 
in frequency bin n, x is the position of the bubble, / is 
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FIG. 2: (Color online) Data from a 553 kg-day WIMP search, 
shown as a distribution in \n{AP) in red. Twenty single nu- 
clear recoil events candidates and 2474 alpha events were ob- 
served. The blue histogram shows the identical analysis for 
data taken in the presence of an AmBe neutron source. We 
define an acoustic cut of 0.7 < AP < 1.3 to select nuclear re- 
coils with an acceptance of 95.8% as determined by the AmBe 
calibration. 



frequency, /mm and fmax are the boundaries of the fre- 
quency band, and pscPj: is the power spectral density for 
the bin with center frequency / for sensor j. The AP 
was scaled to have a value of unity at the peak observed 
in its distribution for nuclear recoils induced by neutron 
sources as shown in Fig. [21 The clear separation seen be- 
tween the alpha peak and the single nuclear recoil peak in 
Fig. [2] illustrates the power of the acoustic discrimination 
to eliminate alpha emitter contamination as a source of 
background for the experiment. 

All data have been subject to a set of data quality 
cuts including the requirement that the chamber expand 
successfully to the desired operating pressure and be sta- 
ble for greater than 30 sec prior to the event. Other 
quality cuts eliminate events with acoustic noise prior to 
the event and events in which the video trigger failed 
to capture the initiation of the bubble. The fiducial vol- 
ume, determined by analyzing the acceptance of the pres- 
sure growth fit cut for events initiated with a neutron 
source, is 92.1 ± 1.8%, equivalent to removing the outer 
2 mm of the liquid volume. This fiducial volume was 
consistent within statistical errors over all neutron cali- 
bration data. The overall efficiency for all data quality 
and fiducial volume cuts is 82.5 ± 1.9%, independent of 
operating temperature. The nuclear recoil acceptance of 
the AP cut alone, shown in Fig. [21 was measured to be 
95.8 ±0.5% in the fiducial volume using the full sample of 
neutron calibration events, resulting in a cumulative effi- 
ciency of 79.1 ± 1.9% for observing a nuclear recoil event. 
Although the standard analysis also identifies multiple 
bubble events with a high efficiency, a complete hand- 
scan of all WIMP search data was performed to ensure 



that none were missed. Therefore, the efSciency for iden- 
tifying muhiple bubble events is 100%. 



IV. BUBBLE NUCLEATION THRESHOLD 

The nuclear recoil energy threshold for the experiment 
was calculated using the Seitz "hot-spike" model of bub- 
ble nucleation[l3| and was benchmarked against calibra- 
tion data. The Seitz model is a two-step thermodynamic 
calculation that begins with the critical bubble radius 
beyond which the bubble will spontaneously grow in a 
superheated fluid: 



2cr 

Pb-Pi = — , 



(2) 



where Pb is the pressure inside of the bubble (vapor pres- 
sure of the fluid), Pi is the pressure outside the bubble 
(expansion set point of the chamber) , a is the surface ten- 
sion of the fluid, and Tc is the critical radius. For bubbles 
smaller than the critical radius the pressure due to sur- 
face tension (right-hand side) is larger than the pressure 
differential across the bubble surface (left-hand side), so 
the bubble collapses. The second step is to calculate 
the enthalpy injection needed to create a critically sized 
bubble, which includes a latent heat term and a surface 
energy term: 



4 f Q 

Et = gTrr^pfc {hb - hi) + A^rl U-T^ 



(3) 



Here pb is the density of bubble vapor, hb and hi are the 
speciflc enthalpies of the bubble vapor and superheated 
fluid, and T is the chamber temperature. In the Seitz 
model, an energy deposition of Et in a volume small 
compared to Tc will nucleate a bubble. All Et and Vc 
values quoted in this paper were calculated using NIST 
REFPROP Version 9.0 [TlLwhich includes models for 
the CF3I equation of state and surface tension (iTj . 

A constant ingress of approximately eight ^^^Rn atoms 
into the chamber per day provided a convenient calibra- 
tion benchmark for the Seitz model threshold, for the 
absolute bubble nucleation efficiency for heavy recoiling 
nuclei, and for characterizing the acoustic signature of 
alpha decays. The decay of one ^^^Rn atom in the cham- 
ber results in three observable events with a readily no- 
ticeable pattern of time correlation driven by the 3.1- 
min half-life of ^^®Po. Figure [3] shows the distribution 
of time differences between 1733 consecutive alpha-decay 
events taken over a period of 4 months compared with 
the results of a fit to a simulation of the expected tim- 
ing of the radon decay chain plus an additional random 
component. The data were best flt by a radon fraction 
of 0.95 ± 0.05, consistent with the expectation that the 
alpha event population is strongly dominated by radon 
decays in the chamber, and unambiguously identifying 
the composition of the alpha event population as equal 
proportions of ^^^Rn, ^^^Po, and ^^'^Po, corresponding 
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FIG. 3: (Color online) Distribution of time differences be- 
tween consecutive alpha-decay events. The solid curve is a 
fit to a simulated time difference distribution, including all 
live time effects and acceptance cuts, based on a component 
arising from decay of ^^^Rn and daughters and a second com- 
ponent arising from random alpha decays with no parent- 
daughter time correlations. The best fit is for a radon fraction 
of 0.95 ± 0.05. For comparison, the dashed gray curve shows 
the expected time difference distribution for uncorrelated al- 
pha decays. The dip in rate around a At of 9 min is caused 
operationally by the forced compression of the chamber after 
a maximum expansion time of 500 sec. 



to nuclear recoil energies of 101, 112, and 146 keV, re- 
spectively. Allowing the bubble nucleation efficiency for 
alpha decays (nuclear recoil plus alpha particle) to float 
as a free parameter in the fit to the alpha time-difference 
distribution yielded a measurement of 100%^2% fo^' ii^- 
cleation efficiency of alpha-decay events at a 15.5 kcV 
threshold. 

By varying the pressure of the chamber, a bubble nu- 
cleation plateau curve as a function of Seitz model thresh- 
old was obtained. The upper graph of Fig. |3] shows the 
plateau curve for single alpha events. The superimposed 
curve illustrates the expected onset of sensitivity to ^^"'Po 
recoils at 146 keV, ^i^Po at 112 keV, and ^'^'^Ru at 101 
keV. A small population of events above the nominal 
nucleation threshold was expected due to the additional 
contribution of the alpha particle to the energy available 
for bubble nucleation. The lower graph of Fig. 2] shows a 
comparable plateau curve for pairs of alpha events sepa- 
rated by less than 500 sec, clearly illustrating the much 
narrower onset of sensitivity to the 101 keV ^^^Rn recoils 
selected by the timing cut. 

The rate of alpha pairs is seen in Fig. U to be con- 
stant within statistical error as a function of threshold 
up to the cutoff, verifying that nucleation efficiency for 
222Rn decays is >75% (at 90% C.L.) up to the ^ispo 
recoil threshold given by the Seitz model. This is con- 
sistent with results from PICASSO [ll], which indicate 
that alpha-decay recoils and ^^F recoils in C4F10 turn on 
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FIG. 4: (Color online) The upper graph is the alpha-decay 
plateau curve for single bubble events, showing rate as a 
function of Seitz model bubble nucleation threshold obtained 
by varying the expansion pressure. The superimposed green 
curve shows the anticipated onset of sensitivity for ^^*Po, 
^^*Po, and ^■^^Rn recoils. The lower graph shows the compa- 
rable plateau curve for pairs of alpha-decay events separated 
in time by less than 500 sec. The superimposed green curve 
illustrates the sharper onset of sensitivity expected from the 
101 keV ^^^Rn recoils selected by the timing cut. In addition 
to the low statistics pressure scan data, the high exposure 
WIMP search data are also included on the plots (the points 
at low threshold with small error bars). 



sharply at the corresponding Seitz model thresholds. 

To determine whether the Seitz model can be extended 
to low energy carbon, fluorine, and iodine recoils in CF3I, 
it is useful to construct the dimensionless quantity 



13 = (rtrack/rc) {pi/pbf'^^ 



(4) 



where rtrack is a measure of the track length of the nu- 
clear recoil in question, and rc is the critical bubble radius 
given by Eq. ([2]) with input conditions (temperature and 
pressure) such that Et as given by Eq. ([3|) is equal to the 
energy of the recoil in question. The ratio of the liquid 
densitiy pi to bubble vapor density pb is used to reduce Vc 
from the critical bubble radius to the radius of the liquid 
volume containing the same number of molecules. The 
distribution of rtrack for a given recoil species and energy 
is found through simulations with TRIM, a Monte-Carlo 
program in the SRIM package that follows nuclear recoil 
cascades in matter [l^-l2l| • The output of TRIM contains 
a list of the spatial coordinates of all displaced atoms in 
the recoil cascade, and rtrack is defined as the square 
root of the maximum eigenvalue of the second moment 
tensor for this distribution of points. For each recoil in 
question, 1,000 tracks are simulated to build the rtrack 
distribution. 

The Seitz model is expected to work well when /? < 1. 
The recoils for which the Seitz model has been verified 
include 6 keV ^^F recoils in C4F10 {/3 = 0.88), 101 keV 



2i8Po recoils in C4F10 {(3 = 0.75), and 101 keV ^ispo 
recoils in CF3I (/3 = 1.02), where the /3 values quoted 
are the median of the distribution. The central 50% of 
the distribution for ^i^Po in CF3I spans 0.86 < /3 < 1.21, 
and the distributions for the other recoils have similar 
widths. For 15(8) keV ^'^'^l recoils in CF3I we find a 
median j3 — 0.70(0.61), supporting the use of the Seitz 
model for bubble nucleation by iodine recoils. Gener- 
ically /3 decreases as recoil energy goes down, i.e., the 
Seitz model should become a more accurate description 
of our threshold as that threshold decreases. 

The situation is less clear for ^^F and ^^C recoils in 
CF3I, which at 15(8) keV have median ^ = 2.02(1.47) 
and /3 = 2.71(2.00), respectively. Previous COUPP cal- 
ibration data '5'] f§\ have shown nucleation rates from 
neutron sources at 30° C to be 50-70% lower than pre- 
dicted by Monte Carlo simulations using the Seitz model. 
Extensive neutron calibration data were taken during 
this run using AmBe and ^^^Cf sources located vari- 
ous distances from the active volume and under varied 
thermodynamic conditions. Each neutron source con- 
figuration was simulated using MCNP-PoliMi [23| and 
GEANT4 [23'! independently, to generate recoil energy 
distribution and interaction rates in the active liquid, us- 
ing the Seitz model in the calculation of bubble nucle- 
ation thresholds. In all cases the predicted nucleation 
rates were larger than those observed, confirming the 
previously observed deviation from 100% nucleation ef- 
ficiency. Given the expected applicability of the Seitz 
model to iodine recoils, we can reasonably attribute the 
observed neutron recoil inefficiency to the ^^F and ^^C 
recoils, with their physically larger energy distribution 
profiles. 

To characterize the observed inefficiency, the data were 
compared to two single-parameter, ad hoc models. The 
first, a "fiat" model, consists of a step function centered 
at the threshold determined by the Seitz theory rising 
to an energy-independent nucleation efficiency, 77 < 1. 
The second model is a function of the energy deposi- 
tion Er and Seitz threshold Et whereby the probablity 
P{Er,ET) of nucleating a bubble is 



P{Er,ET) = 1 - exp 



E — Ex 
Ej" 



(5) 



and a is a parameter describing the width of the turn- 
on. This model has been used by both the PICASSO and 
SIMPLE Collaborations with values of a ranging from 1 
to 1OHJ,01. 

Both efficiency models were fit to the rates of single, 
double, triple and quadruple bubble events for each tem- 
perature set point and several combinations of source 
and source location. The free parameters were rjc,F for 
the flat model and ac,F for the model given by Eq. ([5|), 
and the efficiency for iodine recoils was fixed at 1.0 for 
events above threshold. Both models produced accept- 
able fits as determined by the x^-distribution for Pois- 
son statistics, with the best-fit ryc.F = 0.49 ± 0.02 and 
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TABLE I: Predicted rates for background neutron events aris- 
ing from (Q:,n) reactions and spontaneous fission in the de- 
tector materials near the CF3I volume, and for background 
gamma events from the measured ambient gamma flux. Pre- 
dictions are shown for the three different bubble nucleation 
thresholds, based on a flat 49% nucleation efficiency for car- 
bon and fluorine recoils above threshold and 100% efficiency 
for iodine. The sensitivities to gamma interactions are based 
on in-situ measurements with ^^Co and ^'^''Ba calibration 
sources. 



FIG. 5: (Color online) The observed count rates at the three 
thresholds are for one, two, three, and four bubble events 
induced by an AmBe neutron source. The superimposed 
curves represent the MCNP predictions for the bare Seitz 
model (black) compared to the best fit flat and exponential 
bubble nucleation efficiency models, with 77c, f = 0.49 and 
CiC,F ~ 0.15 respectively. The bare Seitz model clearly over- 
predicts the number of observed counts, especially at high 
multiplicities, and these data do not distinguish between the 
flat and exponential efficiency models. 



ctc'.F = 0.15 ± 0.02 (statistical error bars). The com- 
parison of neutron source data to the MCNP predictions 
is shown in Fig. [S]for both nucleation threshold models. 
For reference. Fig. [S] also illustrates the prediction of the 
bare Seitz model, equivalent to the flat model with a car- 
bon and fluorine nucleation efficiency rjc.F = 1-0 or to 
the exponential model with a very large value of ac.F- 

Note that Eq. ([5]) with ac.F = 0.15 provides a much 
slower rise in nucleation efficiency with energy than has 
been observed by PICASSO and SIMPLE and greatly 
decreases the sensitivity of our detector by cutting into 
the low energy portion of the recoil spectrum. However, 
since the data cannot distinguish between these mod- 
els, WIMP-nucleon interactions limits are presented as a 
band with edges defined by the two efficiency models. 



V. BACKGROUNDS 

While efforts have been made to minimize the neutron 
background from sources external to the bubble chamber 
(both cosmogenic neutrons and those generated by spon- 
taneous fission in the surrounding rock) , a non- negligible 
background is produced internally via both the (a,n) re- 
action and spontaneous fission from the ^38^ g^j^^^ "^^"^Th 
decays in the materials surrounding the CF3I volume. 

A variety of materials used in the bubble chamber were 
screened for their content in U, Th and Po-210, the latter 
an alpha emitter abundantly present in lead-containing 
materials such as the PZT acoustic transducers [23. The 



Nucleation 


Total 


Observed (Predicted) 


Threshold 


Exposure 




Event Counts 




(koV) 


(kg-days) 


Nb = 1 


Nt = 2 


Nb = 3 


Ni = 1* 


7.8 ±1.1 


70.6 


6 (1.0) 


(0.3) 


(0.1) 


2 (0.8) 


11.0 ± 1.6 


88.5 


6 (0.8) 


(0.3) 


2 (0.1) 


3 (0.7) 


15.5 ±2.3 


394.0 


8 (3.5) 


1 (1.0) 


(0.3) 


8 (3.0) 



TABLE II: Observed counts and predicted backgrounds for 
each data set. There is a 79.1% efficiency to detect single 
bubble recoils after the all analysis cuts including the acoustic 
cut described above. Multiple bubble events are identified 
with 100% efficiency by hand-scanning the WIMP search data 
with no quality cuts applied. The final column counts single 
bubbles that survive a 530-sec time isolation cut. 



(Q!,n) and spontaneous fission neutron production rate 
and energy spectrum for each material were calculated 
using the S0URCES-4C [26| code supplied with the mea- 
sured ^■^^U, ^'^^Th, and respective daughter isotope con- 
centrations and the total composition of the material 
in question as inputs, assuming natural abundances of 
any {a,n) target isotopes. These neutron spectra were 
then used to describe the sources in MCNP-PoliMi Monte 
Carlo simulations, and a bubble nucleation rate predic- 
tion was generated for each material in the bubble cham- 
ber that could act as an internal source of neutrons. 

Of the materials considered, most are expected to con- 
tribute less than one event per year in total. However, the 
eight "271 PZT piezoelectric transducers epoxied to the ex- 
terior of the bell jar and the borosilicate glass viewports 
were found to contribute a significant background rate to 
the bubble chamber. Both of these materials are particu- 
larly efficient at generating (a,n) and spontaneous fission 
neutrons, because of their relatively high concentration 
of 238^ ^'^^Th and abundance of light nuclei. 

Table U lists the predicted rates of single and multi- 
ple bubble events at the three operating thresholds as- 
suming a bubble nucleation efficiency of 100% on iodine 
and 49% on carbon and fluorine. At each threshold, we 
predict about 0.012 single bubble cts/kg/day in the de- 
tector from the studied sources. The borosilicate view- 
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ports contribute 73% of this rate, the piezoelectric trans- 
ducers contribute another 25%, with the remainder pro- 
duced by a combination of steel, epoxy and other compo- 
nents. These predictions are subject to a systematic un- 
certainty of 25% arising from the uncertainties in mate- 
rials screening, the MCNP propa gati on of neutrons, and 
from the quoted 18% uncertainty [28| in the results from 
S0URCES-4C. 

The efficiency with which gamma interactions nucle- 
ate bubbles in the detector was measured in situ with 
100 /iCi ®°Co and 1 mCi ^•^■^Ba sources placed inside the 
water shield. At 7.8 keV threshold both gamma sources 
produced an excess of single bubble events, correspond- 
ing to bubble nucleation efficiencies for single gamma 
interactions from either source of 1.4 x 10~^. No re- 
sponse above background was observed at the two higher 
thresholds, providing the limits shown in Table HI The 
gamma ray flux seen by the chamber with and without 
gamma sources was measured by replacing the fused sil- 
ica bell jar with a 1.78 kg NaI[Tl] scintillator. Based 
on MCNP simulations of the NaI[Tl] and CF3I targets, 
the measured background flux in the scintillator corre- 
sponds to a rate of gamma interactions in the CF3I of 
3.4 X 10^ cts/kg/day. Taking the nucleation probability 
to be independent of gamma interaction energy, the re- 
sulting gamma backgrounds or limits thereon are shown 
in Table HI The background from gamma interactions is 
^1/3 the neutron background at 7.8 keV threshold and 
negligible at 11.0 and 15.5 keV. The rate of beta de- 
cays in the CF3I is unknown. Taking the worst-case sce- 
nario of an atmospheric abundance of ^*C, the beta-decay 
rate and resulting background would be 3 times that for 
gamma interactions. 



VI. WIMP SEARCH DATA 

WIMP search data were accumulated between Novem- 
ber 6, 2010, and June 17, 2011, corresponding to a total 
exposure of 553.0 kg-days distributed over three different 
bubble nucleation thresholds. The total effective expo- 
sure for single recoil events given the 79.1% detection 
efficiency described above was 437.4 kg-days. Figure [2] 
shows the AP distribution for all data sets combined, 
compared to neutron calibration data. Twenty candidate 
nuclear recoil events and three multiple bubble events 
were observed, compared to a prediction of 5.3 single nu- 
clear recoil events and 2.2 multiple bubble events from 
the backgrounds described in Sec. [V] 

The numbers of counts observed at the three differ- 
ent bubble nucleation thresholds are provided in Table HIl 
along with the predicted numbers of counts from the 
background simulation. The uncertainty on the Seitz 
threshold is calculated by combining our estimated sys- 
tematic uncertainties on the temperature (1°C) and pres- 
sure (0.5 psia). The largest exposure was at a threshold 
of 15.5 ± 2.3 keV with 394.0 total kg-days of live time. 
Including the 79.1% efficiency for detecting single bubble 



recoil events, the effective exposure was 311.4 kg-days, 
yielding 8 single nuclear recoil events compared to a pre- 
diction of 3.5. At this threshold, we observed 1 two- 
bubble event (with 100% detection efficiency) compared 
to a prediction of 1.0. Because of the generous separation 
observed between alpha particles and nuclear recoils in 
Fig. [21 and because some of the events can be accounted 
for as neutron backgrounds, we do not anticipate that 
alpha rejection failure represents a large fraction of the 
observed single recoil candidate events in the 15.5-keV 
sample. If, however, we interpret all of the 8 events at 
the 15.5-keV threshold as alpha discrimination failures, 
then based on 1733 tagged alpha decays we derive a 90% 
C.L. upper limit on the binomial probability of an alpha 
decay registering in the nuclear recoil signal region to be 
< 0.7%. 

Shorter exposures at 7.8 ± 1.1- and 11.0 ± 1.6-kcV 
thresholds yielded 6 single nuclear recoil events each in 
70.6 and 88.5 total kg-days, respectively. Two three- 
bubble events were observed during the 11-keV exposure. 
The observed single recoil rates at lower threshold are 




WIMP Mass (GeV) 



FIG. 6: (Color online) The 90% C.L. limit for this result is 
shown in blue, interpreting all 20 observed single recoil events 
as WIMP candidates with no background subtraction. The 
band represents the systematic uncertainty in the bubble nu- 
cleation efficiency of fluorine recoils (see Sec. IIV[) . A previ- 
ous COUPP result g] is shown for comparison. The direct 
detection limit from the PICASSO experiment is shown in 
cyan [s^], as well as a controversial limit from the SIMPLE 
experiment in dark green [s^, [sj . Limits on neutralino anni- 
hilation in the sun from the IceCube [ssl ]. magenta, and Super 
Kamiokande [sH, black, neutrino observatories are also plot- 
ted. The indirect detection limits from the neutrino observa- 
tions have additional dependence on the branching fractions 
of the annihilation products. Also shown are limits from col- 
lider searches by CDF ^37] and CMS [11]. The two limits from 
CDF take an effective field theory (valid for a heavy media- 
tor) and a modified theory for a 100 GeV mediator. The CMS 
limits use an effective field theory. The gold region indicates 
favored regions in cMSSM [3!^ . 
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FIG. 7: (Color online) COUPP-4kg limits on spin- 
independent WIMP-proton elastic scattering from the data 
presented in this note are shown in blue. A previous COUPP 
result Ql is shown for comparison. Direct detection limits 
from the XENONIO ;4Q] and XENONIOO [4l[ experiments are 
shown in magenta, and the CDMS experiment [Z^l in black. 
The gold region indicates favored regions in cMSSM [3^. 



by three min. No time clustering is observed in the 15.5- 
keV samples. 

Third, a significant fraction of the low threshold events 
are correlated in time with a bubble in the previous ex- 
pansion. A time isolation cut of 530 secfl^ would have 
eliminated all of the high AP events and all of the time 
correlated events in the 7.8-keV data set, leaving only 
two nuclear recoil events. Further, seven of eight nuclear 
recoil or high AP event candidates that would have failed 
a time isolation cut were specifically correlated to prior 
bubbles occurring very near to the water-CFsI interface 
where a faint but visible ring of unknown residue was 
observed on the inner surface of the quartz vessel. A 
time isolation cut would also have removed three of the 
six nuclear recoil candidate events in the 11.0-keV sam- 
ple but would have no effect on the eight events in the 
15.5-keV sample, leaving two, three, and eight nuclear 
recoil candidate events in the 7.8-, 11. 0-, and 15.5-keV 
samples respectively. These numbers of counts are still 
higher than the 0.7, 0.8, and 3.5 events predicted by our 
neutron simulation, but the significance of the excess is 
diminished by the lack of any method for estimating the 
fraction of the spurious events which still pass a time 
isolation cut. 



significantly higher than the 0.7(0.8) events predicted by 
the neutron simulations at the 7.8(11.0)-keV thresholds, 
suggesting an excess of single nuclear recoil events in the 
7-15 keV range. 

We note however that this low threshold population of 
candidate nuclear recoil events differs in three ways from 
what would be expected from true single nuclear recoils. 
First, the AP distribution for the single nuclear recoil 
events in the low threshold samples is noticeably broader 
than was observed in calibration neutron events taken 
under the same operating conditions and has a signifi- 
cant tail to higher values of AP. This can been seen in 
Fig. [5] Whereas the nominal AP cut has been measured 
to be 96% efficient for calibration neutron events, relax- 
ing our AP cut to 0.7-1.5 increases the number of nuclear 
recoil candidates from 6(6) to 10(8) in the 7.8(11.0)-keV 
samples. The AP distribution for the 15.5-keV sample is 
consistent with the neutron calibration data. 

Second, a significant fraction of the events in the 7.8- 
keV sample occur in statistically unlikely clusters. Using 
the less restrictive 0.7-1.5 AP cut, and additionally con- 
sidering events with acceptable AP but narrowly rejected 
for other data quality cuts, we obtain a sample of 12 nu- 
clear recoil candidate events or near misses distributed 
over a period of 14 days. Three of the 12 events occur 
in a 3-hour time period, with two occurring eight min 
apart. A second group of five events occur in an 8-hour 
time period, with three events occurring in a ten-min in- 
terval. Two events in the 11.0-keV sample are separated 



VII. CONCLUSION 

Because a time isolation cut was not benchniarked 
prior to our low background running and given the sys- 
tematic uncertainties in the neutron background simula- 
tions, no background subtraction has been attempted. 
Our limits are therefore based on treating all 20 nu- 
clear recoil events passing our cuts as dark matter can- 
didates. The resulting 90% C.L. limit plots for spin- 
dependent WIMP-proton and spin-independent WIMP- 
nucleon cross sections are presented in Figs H] and [71 re- 
spectively. The calculations assume the standard halo 
parameterization [3^], with = 0.3 GeV cm~'^, 
Vesc = 544 km/s, ve = 244 km/s, vq = 230 km/s, and 
the spin-dependent parameters from the compilation in 
Tovey et al.[31]. 
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